Random walks of molecular motors arising from diffusional encounters with 

immobilized filaments 
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Movements of molecular motors on cytoskeletal filaments are described by directed walks on 
a line. Detachment from this line is allowed to occur with a small probability. Motion in the 
surrounding fluid is described by symmetric random walks. Effects of detachment and reattachment 
are calculated by an analytical solution of the master equation in two and three dimensions. Results 
are obtained for the fraction of bound motors, their average velocity and displacement. The diffusion 
coefficient parallel to the filament becomes anomalously large since detachment and subsequent 
reattachment, in the presence of directed motion of the bound motors, leads to a broadening of the 
density distribution. 

The occurrence of protofilaments on a microtubule is modeled by internal states of the binding 
sites. After a transient time all protofilaments become equally populated. 
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I. INTRODUCTION 

Molecular motors are proteins that convert the free en- 
ergy released from chemical reactions into directed move- 
ments [2, 2] ■ Here we focus on linear cytoskeletal motors 
which move along cytoskeletal filaments. Many of these 
motors are processive in the sense that a single motor 
molecule can move a cargo over a large distance. The 
most prominent examples are (conventional) kinesin and 
certain types of myosins, which move along microtubules 
and actin filaments, respectively. In the cell, these mo- 
tors are involved in transport processes, reorganization 
of the cytoskeleton and cell division _1J. However, exper- 
iments on the movements of molecular motors can also 
be done in vitro, which has lead to the development of 
various single molecule assays. In these experiments, one 
can measure the velocities, step sizes, walking distances 
and forces for single motor molecules, see, e.g. 0- In 
additions, they have stimulated a lot of theoretical work 
devoted to the walks of molecular motors along filaments, 
see, e.g. |3- 

In the experiments, one observes that even processive 
motors unbind from their filamentous tracks after a cer- 
tain walking distance, which is typically of the order of 
1 ixm. For a kinesin molecule, this means that it makes 
about 100 steps of 8 nm before unbinding nil. Myosin V 
motors have a comparable walking distance, but a larger 
step size of 36 nm, so that they detach after about 30-50 
steps [5 ^7^. Unbound motors then diffuse in the sur- 
rounding fluid until they rebind to the same or another 
filament and continue their directed walk. 
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On larger scales, the motors thus perform complex ran- 
dom walks, which consist of alternating sequences of di- 
rected movements along filaments and non-directed dif- 
fusion in the surrounding fluid, as shown schematically 
in Fig. ^ These random walks have been discussed by 
Ajdari using scaling arguments Recently, we have in- 
troduced lattice models to describe the random walks of 
the motors as random walks on a lattice, where certain 
lines of lattice sites represent the filaments 0, ^3 ■ When 
bound to these lines, the motors perform directed ran- 
dom walks. Detachment from these lines is allowed to 
occur with a small, but non-zero probability. Diffusive 
motion in the surrounding fluid is described by symmet- 
ric random walks. 

These models are designed to study generic properties 
of motor movements, but they can also be used to de- 
scribe specific motor molecules, since all model param- 
eters can be determined from the measured transport 
properties, see Ref. @- In addition, motor-motor in- 
teractions can be easily included into these models, e.g., 
mutual exclusion of motors from the binding sites of the 
filaments, which leads to self-organized density profiles 
in closed systems |^ and boundary-induced phase tran- 
sitions in open tube systems {lH. 

For the random walks of single motors or, equiva- 
lently, for an ensemble of non- interacting motors, we have 
obtained a number of exact results for the cases of a 
single filament embedded in two-dimensional or three- 
dimensional space as reported in jlflj. In particular, these 
random walks exhibit anomalous drift behavior; the av- 
erage position of the motor advances slower than linearly 
with time. The same drift behavior is found for the move- 
ment along a single filament immobilized in open com- 
partments with the same dimensionality, which are more 
easily accessible to experiments j^] . In the present article, 
we present a detailed derivation of the analytical results 
of Ref. jlOl | for movements in two and three dimensions 
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without confining boundaries. 

Analytical results are obtained by the following 
method, which is a variant of the method of Fourier- 
Laplace transforms for random walks in homogeneous 
space, see, e.g., 12, 13J: By using Fourier-Laplace trans- 
forms of the probability distributions, the master equa- 
tions of the random walk can be transformed into a set of 
algebraic equations, one of which, however, requires the 
evaluation of a non-trivial integral. Solving these alge- 
braic equations, solutions for the Fourier-Laplace trans- 
formed probability distributions and their moments are 
obtained and closed expressions in terms of integrals can 
be derived for the time-dependent probability distribu- 
tions and moments. These can, on the one hand, be eval- 
uated numerically to obtain results for all times; on the 
other hand, asymptotic results for small and large times 
can be obtained fully analytically by using the Tauberian 
theorems, which we summarize in the appendix. In this 
way, we derive expressions for the fraction of bound mo- 
tors, the average displacement and dispersion, and the 
effective velocities and diffusion coefficients. The ana- 
lytical results are compared to data from Monte Carlo 
simulation and are found to be in very good agreement. 

In addition to the anomalous drift behavior, the ran- 
dom walks of molecular motors also exhibit strongly en- 
hanced diffusion in the direction parallel to the filament. 

Our article is organized as follows: We start with the 
two-dimensional case in section Inland discuss the three- 
dimensional case in section imi In both cases, we derive 
probability distributions and their moments for both the 
bound and unbound motor. The fact that filaments may 
consist of several protofilaments is taken into account in 
the final subsections of sections ^ and IIIII where these 
are modeled by several internal states of the bound mo- 



P„,m(t+1) = 

Pn,o{t+l) = 

PnAt+l) = 

P„,_l(t+1) = 

As initial condition we take an ensemble of particles at 
n = m = 0, so 



Pn,m{t = 0) = 5nfiSm,0- (6) 



Let us now define the Fourier-Laplace transforms of the 
probability distribution along the filament Pi,(n,t) = 



tors. In section HVI we extend the discussion to include 
a variable sticking probability for motors arriving at the 
filament. At the end, we include a short summary of our 
results. 



II. RANDOM WALKS IN TWO DIMENSIONS 

Consider a discrete time random walk on a two dimen- 
sional square lattice with lattice sites labeled by integer 
coordinates (n,m). At each step a particle has proba- 
bility 1/4 to jump into any of the four directions. For 
modeling the motion of a motor on a microtubule, we 
choose different behavior on the line with m = 0. Here 
the probability to jump from {n,0) to (n,±l) equals je, 
while jumps to {n + 1,0) have probability 1 — 7 — -^(5 — ^e; 
the probability to jump to (n — 1, 0) is i(5 and the prob- 
ability to make no jump is 7, see Fig. |2l The latter pa- 
rameter is needed for modeling realistic situations, where 
the diffusion coefficient in the fiuid is much larger than 
on the filament 0- The ordinary random walk in two 
dimensions has 7 = 0, S = 1/2, and e ~ 1. We shall 
assume that the escape probability e is small. For e = 
the problem amounts to a directed random walk on the 
line with m = 0. The average speed of a motor particle 
on the filament line isfb = 1 — 7 — (5— 5 e. Per step there 
is a probability ^ e to unbind. The probability that the 
motor is still bound after t steps is 

(l_i,)*«exp(-iei). (1) 

The master equation for this dynamics reads 

(2) 
(3) 
(4) 
(5) 

I 

Pn.o{t) and the full distribution P„.m(i), 

and 

^ ^ piqva-\-irn 
t—{) m,n— — 00 



I 

^ Pn+l,m + ^ Pn-l,m + ^ Pn,m+1 + ^ Pn,m-l [m 7^ 0, ±1) 

\ PnA + \ P«,-i + (l-7-^e-^ <5)P„-i,o + + 7Pn,0 

^ Pn+ia + ^ Pn-1,1 + ^ Pn,2 + 

1 1 1 e 

^ + ^ + ^ Pn-2 + ^Pn,0- 
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Then the master equations are reduced to an algebraic equation relating Pi^(r,s) and P{q,r,s) as given by 
I 



(1 + s — - cos '^'^^ ''")P{Qj s) = 1 



+ 7- 



1 -e 



■ COSQ 



Phir, s). 



(9) 



Here the first line is what one would get in the case of a symmetric random walk in two dimensions, and the second 
line corrects those terms that are changed by the presence of the filament. This equation has the obvious solution 



P{q,r,s) = 



1 + [7(1 — cosr) + i (1 — e)(cosr — cosg) + iv]^ sin r]Pi,{r,s) 



s + 1 — ^ cos q 



cosr 



(10) 



By integrating this result over q we also obtain Ph{r, s) on the left hand side. It thus satisfies a linear equation, that 
can be easily solved. Introducing the variable fi via 



we may use the equalities 

'■^^ dq 



cosh /i = 2 + 2s — cos r or sinh ^ = \/ (2 + 2s — cosr)^ — 1, 

dq 



1 



1 



and 



cosq 



/q 27r cosh n — cos q sinh /i Jq 27r cosh fi — cos q sinh fj, 

After some computation, we then end up with the probability distribution 



(11) 
(12) 



s + 1 - 7 - (1 - 7 - ^ - i 



1 



s + (1 — 7)(1 — cosr) + 56 (cosr — e ^) — iv^, sinr 
I 



(13) 



for the motors bound to the filaments. The probability value at r = gives us the Laplace transform No{s) of the 

distribution for all, bound and unbound motors follows probability A^o(0 = X]n^n,o(i) that the motor particle 

via Eq. (|10|) . It is easy to check that this is correct for is bound to the filament line with to = 0: 
e = (random walk in one dimension) and for 7 = 0, 
(5 = i , and e = 1 (non-biased random walk in 2d). 



A. Properties of the motors bound to the line 

1. Survival fraction 

In the following, we extract the transport properties of 
the motor's random walks from the solution H13|) . The 



The inverse is 



Noit) 



-.^l + sY+i 5+i£(i_e-P) (1_ e)s+e^s{l + s) 



. ^ ^ ds ,t , , dx e(l-a;)*+i/2 

iVo(0=c^^(i + .)*iVo(.) = y^ VT^IT^TTTT^ 

I 



dy (1- £^ty)*+i/2 
Vvy[l + (l-2e)y]- 



(14) 



(15) 



This expression is exact. It holds for all t and for all obtained in this way are shown Fig. |31 (lines) for three 
e, S and can be evaluated numerically. Values for No{t) 
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different values of the detachment rate e. Comparison to 
results of Monte Carlo simulations (data points) shows 
that the agreement is very good. 

Let us now derive the asymptotic behavior for small 
and large times, respectively. If e is small and i » 1, we 
have 



Nois) 



1 



s + 



(16) 



Then the inverse behaves as 



Noit) 



ds 



27ri s - 



dy 



(17) 



The second integral is an obvious limit of the last integral 
in Eq. (|15ll . For deriving the short-time behavior we may 
start from the series of equalities as given by 



-yeH 



1 * - 1 

+ - 



y+l 



1 



y + l 



y 



e^t e-y""* + ye^t - 1 
y + l 2/2 



+ 



(18) 



The integrals over the exponential terms are most easily carried out using dimensional regularization. To show how 
that works, let us consider the last term. We need to consider the expression 



1 



\yy^-^'^{e-y'"' + yeh-l) = -^' 



2i\3/2-r 



r(n- -) 



{n-l){n-\) 



r 



(19) 



in the limit of small n, where we used that, in dimensional 
regularization, integrals of powers are put equal to zero. 
The limit n — > can now be taken. Using the same 
procedure for the other terms, we obtain 



I 



+ eH 



30F 



(20) 



for small t, which represents a series in powers of e\/i. For 
t <C this is somewhat surprising: although the mo- 
tors detach at times 1/ e ^ the recurrent behavior of the 
random walk brings them mostly back to the filament. 
For t ^ this just says that the motor did not have 
enough time to escape from the line. The half-integer 
powers are related to the long range diffusion away from 
the filament. 

Note that, for the short time limit, we have assumed 
that t ^ I, i.e., that t is large compared to the time 
required for one step of the random walk. In that interval 
the approximation H16|l holds with the term e y/s being 
due to diffusion. 

For smaller times, the random walk exhibits discrete 
steps. Our short time result would hold for arbitrarily 
small times in the limit in which the walk on the line 
becomes a continuous time random walk. A continuous 
time random walk with exponential waiting time distri- 
bution is obtained approximately, if 1 — 7 is small, which 
is the case for realistic applications of our model to molec- 
ular motor setups |9| . The same remark holds for all the 
short time results discussed in the following. 

Another way to derive H2U|I is to expand (|14|) in powers 
of l/^/s and to use the Tauberian theorem which states 
that the inverse Laplace transform of No{s) = as~°' is 
given by 



T{a) 



a-l 



(21) 



see appendix IXI and, e.g., Ref. This theorem holds 
both for positive and negative values of a. It also shows 
that positive integer powers of s of Nq{s) in the limit of 
small s do not contribute to long time tails. 

For large times, the expression lfT7|l for No{t) 

can be evaluated using the expansion l/{l+y) w l—y+y^ 
for small y which leads to 



No{t) 



I 



(22) 



Alternatively, one can expand (|14() in powers of ^/s and 
use again the Tauberian theorem. It follows from H22|) 
that for large t the probability to be bound to the filament 
decays as i^^/^ jj^ agreement with scaling arguments 



2. Average position and speed on the filament line 

The expression (|13f) for the Fourier-Laplace trans- 
formed probability distribution P\^{r,s) contains much 
more information. At linear order in r we get the Laplace 
transform of the average position of motor particles along 
the filament line, 



d 



r=0 



We obtain 

N,is) = v^N§is) 



[(1 - e)s + eyi(lT^]2 
In the limit of small e and large t, this implies 



(23) 



(24) 



(25) 
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We invert the Laplace transform by taking the 1/s term 
apart. We then get for the average position the asymp- 
totic behavior 



Ni{t) 



_ 2^^^ p _dx_ e-"* 



e 



dy 1 



VV (1 + 2/)' 



(26) 



for large t. The exact expression is 



Niit) 



dy 



Vy[l + 2/(l-2e)]2 



(27) 

which deviates from 1)26(1 for times of order unity. The 
full expression H27|l is evaluated numerically and plotted 
in Fig. ^ The same figure contains data points as ob- 
tained from Monte Carlo simulations which confirm the 
analytical result (|27|l . 

For short times we proceed as above. We expand 



Niit) 



2vb 



dx 



xeH 



3 x/TT 



(28) 



This leads to A'^i « v\jt. Thus, the average position fib 
and speed Vb of the motors bound to the filament are 
given by 

Ni(t) 2 eVt, 

No{t) 3 ^/TT 



and 



Vh 



drib 
dt 



^ : 



(29) 



where is the average speed if the particles did not leave 
the line. 

For large times, the asymptotic expression (|26|l leads 

to 

m{t) « ^(i-4t) 



and lib(t) « !^!i^(l-^). (30) 
e ey/irt 

Therefore, for large t, the average speed of bound motors 
behaves as 



(31) 



The last relation confirms the scaling Vh{t) ~ VhNo{t), 
which has been used in the scaling approach |^, ^ . The 
effective motor velocity is reduced by a factor ~ No{t), 
i.e., by the probability that a motor is in the bound state. 
The relation vy, NqVi, also applies to a simple two- 
state random walk, where motion is directed in one of 
the states only. In contrast to the simple two-state ran- 
dom walk, however, the probability Nq is time-dependent 
here. The factor tt/2 in is solely due to the fact that 
only the bound motors are considered. We will show be- 
low, that this factor is absent, if all, bound and unbound, 
motors are considered. 



3. Dispersion and diffusion coefficient on the filament line 

By expanding the expression H13|l for the Fourier- 
Laplace transformed probability Ph{r., s) up to second or- 
der in r, the second moment of the distribution of bound 
motors is found to behave as 



4Vs + s2 



)N^{s). (32) 



The Laplace transform is again inverted by complex inte- 
gration. In the contour integral for N2{t), we may replace 
e** by e"* — 1, since the subtracted integral vanishes, as 
can be seen by closing the contour of integration in the 
right half plane. 

Closing the contour along the negative real axis, we 
obtain 



2(1-7- e) 



da; 



[(l-a:£')*-l]Vl- e^x 2vl 



n € ■ 



dx 



[l + (l-2e)a;]3V^ 
[l-(l-a;e2)t]Vl- e^a 



da; 



[(l-a;e^)*-l](l- e^xf/^ 
[l-F(l-2e)a;]3a;3/2 



An e • 



dx 







[l + (l-2e)a;]2V^ 

(1 - 2 e 2x) [(1 - xe^)* - 1 -f e ^xt\ [1 - (1 - 2 e )x\ 
[1 + (1 - 2 e )a;]2a;3/Vl - e^a; 



2z^ (1-e— *)(l-3a;) 2(1-7) 

7re4 jo (1 -t- a;)3a;3/2 



da; 1 



e2 7o TT (l + a;)27i 
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1 

4^2 



dx (1 — e 

TT 



(1 + x)^x^/x 



(33) 



We have evaluated this exact expression numerically and compared it with simulation data. As shown in Fig. (HJ the 
agreement is again very good. 

For short times we proceed as above. 



N2it) 



2vl 



dx 



{xe't~\x^eH''){l-?,x) 
(l + a;)3x3/2 



2 2/, 87-7e^^/^ 

V (1 - "T —r) + (1 - 7)^ + — ^ ^• 

5 VTT 6 ^-K 



1 + xe^i 



,7/2 



(34) 



r 



This implies the normalized second moment 

n\ = ^rr^TV ~ <tHl - - ^) + (1 - j)t + ^ ^ 

(35) 

and the dispersion 



Anl 



15 ^TT 3 ^TT 

(36) 

This result holds again in the limit of continuous time, i.e. 
for 7 close to one. The leading term taking into account 
the time discretization is (1 — ^ — V\,)t, i.e. the dispersion 
at short times, which arises from the walks along the 
filament, is smaller in discrete than in continuous time. 
For the large time result, which we will derive next, the 
choice of continuous or discrete time makes no difference. 
The relative dispersion 



^2 



2 

15 



(37) 



is small since t <^ l/e^- Using the dispersion An^, we 
may also calculate the time-dependent diffusion coeffi- 
cient 



^(1-7) + 



^ ,3.3/2 



2-1-1 ey/t 



4 ^ 
(38) 

In the scaling regime t ^ Dt,{t) is much larger than 
its limiting value i (1 — 7) for e = 0. This enhanced 
diffusion arises from the fact that each motor may detach 
from the microtubule with probability exp(— r) for any 
value of r = ^ e according to Eq. (Q). This leads to a 
considerable broadening of the bound motor distribution. 

For large t one can make a change of variables in the 
integral expression H33|l and use y = e ^Tx as the new in- 
tegration variable. This leads to the asymptotic behavior 



N2{t) 



M^a44, + l^a-^) ,38, 



e^^' 2 eVt' 



e\/7rt 



and 



Anl^4i^-n-'^)t+ (40) 



The diffusion coefficient behaves as 

^^.,(l-7)^/^ 



V 

2e 



\(4 



TT - 



4e^/^ 



(41) 



for large t. The limiting value of the diffusion coefhcient, 
-Db(oo) ~ (wb/e)^ is large compared to the diffusion co- 
efficient of the one-dimensional random walk along the 
filament, -Db(O) = (1 — 7)/2. This broadening of the dis- 
tribution occurs since the unbound motors lag behind the 
bound ones, which implies that also the rebinding motors 
lag behind those that have been bound for some time. 



^. The density profile on the filament 

At large times n'^h scales as and one may look for 
a scaling form of the density. For small s and small r, 
the Fourier-Laplace transformed probability distribution 
Ph{r,s) as in (|13l behaves as 



Phir,s) 



(42) 



The inverse Fourier transform follows from the contin- 
uum limit of Eq. (jT)). It then follows that, in this limit, 
Pn,o{s) = for n < 0, while for positive n, one finds 



_Lp-"(«+«Vs)/»'b 



Vh 



(43) 



the inverse Laplace transform of this expression now leads 
to Pn,o{t) = for n > v^t which implies that the overall 
motion is slower than ballistic. For n < v\,t, we obtain 



1 

Pnflit) « / ds 



(44) 



After a change of variables from s to w = ^/s, the it- 
integral can be taken over the whole real axis (at the 
expense of a factor i ) . Evaluation of this Gaussian in- 
tegral leads to 



Pnfi{t) 



exp 



4ub(wbi - n)^ 



(45) 
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for n > 0. This expression vanishes hnearly as '--^ n for 
small n and even exponentially fast as n approaches vi,t 
from below. The density profile Pn,o{t) as given by lH^ 
is plotted in Fig.[Slfor several values of the time t. Com- 
parison with the results of MC simulations shows that 
the asymptotic expression as given by H45|l is very good 
for times which excess about 8000 time steps. At smaller 
times, the asymptotic expression overestimates the maxi- 
mum of Pn,o and underestimates the tails of Pn,o for large 
n, see Fig.|Sl Simulation data are obtained by averaging 
over 5 x 10^ realizations of the random walk. 

It is somewhat tedious to show that the moments 
A^'o, iVi and N2 as obtained from H45I) agree with the 
previously derived expressions. To verify this, one 
may, in an intermediate step, use the substitution n = 
2vht{\/ y + 11^ — y) which implies 



Pnfi{t) An 



ty 



e^tdy 



Ax 



■'(x+y+xy)e'^t 



(46) 



B. Properties of the unbound motors 

Eventually every motor will unbind and diffuse in the 
surrounding fluid. We now discuss the effects of the fil- 
ament on the behavior of the unbound motors. The 
Fourier-Laplace transform P{q, r, s) of the probability 
distribution Pn.m{t) of bound and unbound motors as 
given by Eq. H1U|) can be rewritten as 



P{q,r,s) = Ph{r,s) +Pnh{q,r, s). 



(47) 



The first part is Ph{r,s) as given by Eq. H42I) . describ- 
ing the bound motors. The second part describes the 
corresponding probability distribution 



Pub(<Z,r,s) = 5] ^e'*™+'"-"P„,„.(s) 



s + 1 — ^ COS q — h COS r 



7(1 — COS r) + (cos r — cos q) + iv\, sin r 



s + 1 — \ cos q — J cos r 



Pb(r,s). 



(48) 



of the unbound motors. For small s, small r, and small q, taking into account that r ^ s and q ^ ^/s, the expression 
(|?7jl leads to 



P{q,r,s) = Pb {r,s) + Pub iq,r,s) 



1 



s — ivt,r + e-y/s (s — iv^r + ey/s){q'^ + 4s) 
I 



(49) 



Position and longitudinal diffusion 



which leads to 



Expanding H48() or, for large times, (|49|l in powers of 
r and q yields the moments of the distribution of un- 
bound motors. The total fraction of unbound motors 
as obtained from (|47|1 for r = g = is, of course, 
NQ^{t) = I — No{t) expressing motor conservation. The 
average longitudinal position of the unbound motors is 
given by 



N^^'is) EE ^^nP„,™(s) = -z|-Pub(g,r,s) 

n m^O 



VhNois) 



Vi(s + e^/^)2 



r=q=0 

(50) 



as follows from (|48|l and (|49|l . The last equality holds for 
small e and small s, i.e. for large t. 

For small t, one has from the exact expression in (|50|l 



3/2 



30F 



.(1_ eVTTt), 



(51) 



and 



n^h ~ 2«bi(l - ^ eVTTt) 

2 3 
Vuh ~ x«b(l - -T e\/7rt) 
o o 



(52) 



for the average position and speed of the unbound mo- 
tors. For large t, the asymptotic expression given in H50|) 
leads to 



TTe^ 



00 

As 



which implies 



and 



jT-ub 



3/2(s+l)' 



0re 2 e^/t 

2ub Vt 



(53) 



(54) 
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Whereas each individual motor has zero average velocity 
in the fluid, the statistical velocity is non-zero, since 
it is driven by unbinding from the cloud of motors moving 
on the filament: The cloud of unbound motors advances, 
because motors rebind to the filament and others detach 
from it, and those detaching have propagated a certain 
distance compared to those rebinding. 

Since for large times all motors are detached from the 
filament most of the time, equation (|54|l gives the asymp- 
totic displacement of a motor, if averages are taken over 
all bound and unbound motors, i.e., 



The last relation implies that, in the longitudinal direc- 
tion, the released motors spread in a broad cloud, not 
narrowly centered around its average. 
At long times, one obtains 



Vh 



(55) 



at large times as follows from (PU, and N^^ « 
1. The effective (time-dependent) velocity at large times 
is v{t) « VhNo{t), as predicted by scaling arguments 8, 
1^ . This can also be obtained by inspection of Eq. (|Sn|l : 
The displacement of all motors is given by 



^i(s)-fiVi"b(s)=t,b 



(56) 



The normalization factor for all motors is unity. The 
effective velocity is obtained from the inverse: 



vit) 



Vh- 



dt / 2TTi 



ds 



= ^'b / —e''No{s)^VbNo{t) 
(57) 

Fig. 21 shows the displacement of the bound and un- 
bound motors, obtained by numerical evaluation of the 
exact integrals (dotted lines) and corresponding simula- 
tion data (circles for the bound motors and diamonds for 
the unbound ones). The displacement averaged over all 
motors, n = UhNo+fiuhNQ^ , is also shown in Fig.0](sohd 
line and squares). For small times it is equal to the dis- 
placement of the bound motors, for intermediate times it 
interpolates to the curve for the unbound motors. Since 
almost all motors are detached for large times, the dis- 
placement of all motors is then equal to the displacement 
of the unbound ones. 

The longitudinal position has the second moment 



Zs^ s 



Vs(s -t- eys)3 ■ 
At short times, this means 



iV2"^i) 



15V^ 



16 



(58) 
(59) 



and thus 



and An^t) « vlf±il-l'-2^). 



(60) 



ub 



and 



(61) 



This implies the dispersion and the diffusion coefficient 
as given by 



2{tt - 2)vlt 



ub 



and 



2)v- 



(62) 



The order of magnitude of the diffusion coefficient is 
again w^/e^. The prefactor 1 — 2/7r~0.36is slightly 
smaller than the prefactor 2 — 7r/2 ~ 0.43 of Eq. ||1T|| . 
The order of magnitude -Dy ~ u^/ tells us that, like 
on the line, longitudinal diffusion is strongly enhanced 
by the unbinding from and binding back to the line. 



2. Transverse diffusion 

The diffusion behavior perpendicular to the filament 
is determined by H49|l up to quadratic order in q. The 
average transverse position vanishes and the dispersion 
of the transverse position is given by 

A..^(,s) = ^is) = ^ - i^iVo(.) « ,^3/2(;+,^) 



and 

Am^{t) =m?(t) 
This implies 



-St 



St 



Am^{t) 



2e^3/2 
30F 







and 



s3/2(s + £2) 



20F 



(63) 
(64) 

(65) 



at short times. It is small because the motors started in 
a state bound to the line. At large times one has 



1 



(66) 



The transverse diffusion coefficient 

D^{t)^]-{l~^) = ^[l~No{t)] (67) 
4 eVirt 4 

approaches the free value j. So the transverse diffusion 
starts out very small, and finally reaches its value in free 
space. The quantity Am^(i) is plotted in Fig. [7| There 
are small deviations in comparison to the simulations at 
small times arising from the time discretization. 
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3. Spatio-temporal density profile of the unbound motors 

Finally, we derive the density profile of the unbound 
motors. After inversion of the Fourier-Laplace trans- 
forms, the profile H49|l becomes in real space 




Comparing with H43|) we have a shifted value for e n and 
an extra factor e . The temporal form thus becomes 

e{en + 2\m\vb) ( {en + 2\m\v\,f \ 

2^{v^t - n)3/2 1, Avy^iv^t - n) ) 

(69) 

for m 7^ and n > 0. This again vanishes for n —> Vbt. 
For n it is finite when in is finite, while it is zero 
for negative n, in the present approximation. Though 
detachment followed by diffusion to negative n may oc- 
cur, the attainable values are of order unity, and vanish 
to leading order in Wb/ e^- 

C. Pinning line with several tracks 

Microtubules, the filament tracks of kinesin and dynein 
motors, consist of 13 parallel protofilaments 0, each of 



which provides a possible track for these motors. To 
incorporate this in our model, we assume that the pinning 
line at TO = has k internal states where k may be equal 
to 13. The average occupation of each of these states is 
denoted by p;^ q with j — 0, .., — 1. There is a small 
probability, ^ C, that a motor goes from track j to j + 1 
within one time step, and similarly for going to track 
j — l. To take into account the cylindrical structure of the 
microtubule, we identify j = k with j = 0. We assume 
that after detaching, the motor may randomly go to the 
right or to the left of the tubule; likewise, when attaching 
to the tubule either from right or left, the motor has equal 
probability to attach to any of the tracks. 

In this model the total fraction of motors at position 
n along the tubule is 

k 
J = l 

The dynamics given by Eqs. ® > Oi and Q remains 
valid and leads to the same solution, as it is insensitive 
to the internal distribution over the tracks. 

The motion on the individual tracks is described by a 
master equation analogous to Eq. (j2Jl and given by 



+ {{P^V+Pn:o') + {l~OPl,- (71) 

When summed over j this indeed leads back to Eq. ©. To solve Eq. H71|l . another Fourier transform is needed which 
is defined by 

k 

P^o^T^Ploe''^ (72) 
i=i 

with u! = 2TT£/k and £ — 0, .., k — 1. We assume that, at t = 0, all motors are located on the track j = which 
corresponds to the initial condition 

Plmit = 0) ^ Sn^oSr,.,oSjM- (73) 

Going to the Fourier-Laplace transform, wc find 

2 

" ^ ' ^ 2 + 2s- (2-27-(5-e)e'"--27-5e-"^ + 2C(l-cosw)' ^ ' 

The surviving fractions on the individual tracks are found from the expressions for r = 0. In the time representation, 
the surviving fractions behave as 

pi{t) ^Y^PnAt) = T^oit) + I e-^^'"e-5^*-^(^-™^")*. (75) 



The first term describes the symmetric distribution of the motors over the k tracks; in Eq. H17|) we showed that it de- 
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cays algebraically. The second term represents the asym- 
metry arising from the initial condition that only one 
track is populated. There are two decay mechanisms of 
this asymmetry. The term exp(— ^ et) expresses that de- 
tachment and subsequent reattachment to randomly cho- 
sen tracks restores the symmetry. This happens in partic- 
ular for motors that quickly return to a randomly chosen 
track, implying ordinary exponential decay. The factor 
exp(— (^(1 — cosuj)t) expresses that hopping to neighbor- 
ing tracks also restores the symmetry. For large k, the 
smallest of these rates is 27r^(^/A:^. So the exchange pro- 
cess dominates when this value is larger than e/2. 

Therefore, the asymmetry between the average occu- 
pation of the various tracks disappears after the discussed 
transient time. 



mi = TO2 = represents the filament which attracts 
and binds the motors. Away from the filament the jump 
probabilities are equal to ^ , while on the filament they 



and 2 S in the forward and 



are given by 1 — 7 - 
backward direction, respectively, equal to ^ e for each 
of the four sideward directions, and equal to 7 to make 
no step, see Fig. |51 The average short-time velocity is 
— I e while the sticking probability is 



exp[ 



7 



t]. 



III. RANDOM WALKS IN THREE 
DIMENSIONS 



Now, let us consider the same kind of random walk 
on a three-dimensional cubic lattice, in which the line 



We denote the transverse coordinate as m = (mi, 1122) ■ 
The master equation for this dynamics has the form 



J 



P„,o(i+l) 
-Pn.p(i+1) 



(m^O,p) 



1 2 11 

p 



(76) 

(77) 
(78) 



P'(#P) 



r 



In these equations, p and p' denote the four transverse three vectors. We can follow the same steps as in the two- 
nearest neighbor vectors which connect a filament site to dimensional case using again Fourier-Laplace transforms. 



the four adjacent non-filament sites, p — (0, ±1), (±1, 0) 
The summations over p in l|76|) and H77|l or p' run over 
the four possible values. Eq. H78() holds for any of the four The equivalent of Eq. (|10() becomes 
values of p with the sum over p' running over the other 



The Fourier transformation in the perpendicular direc- 
tions leads to a transverse Fourier vector q = (91,92)- 



P(q, r, s) 



3 + [37 + i (5 - 67 - 4e - 3(5)e*'- - i (1 - 3(5) e"" - (1 - e)(cosgi + cos 92)] Ph{r, s) 



3 -I- 3s — cos r — cos qi — cos 92 
By doing the integrals over qi and q2, we derive the expression for Ph{r, s). The necessary integral is 

1 



(79) 



/(r, s) 



'"dqi /•2-dg2 



2tt Jq 27r 3 -|- 3s — cos r — cos qi — cos 92 
"'x(m) 



with m : 



(3 + 3s — cosr)2 ' 



(80) 
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where K{m) — dif)/\/l — msin^ is the complete elliptic integral of the first kind. We also use the relation 
'-dq, r^dq, COS cos 52 (3 + 3. - COS .) - 1, ( 



/q 27r Jq 27r 3 + 3s — cos r — cos qi ~ cos q2 
where /(r, s) is given by (jHOJ. Wo then get 

P^(r s) = (82) 

e+[3(l-e)s+i(e-3(5)(e-"--l)-i(6-67-3(5-5e)(e''--l)]/(r,s)' ^ ' 



r 



For large s, one may use I{s) ~ l/(3s), to verify that 
Ph{r, s) ~ 1/s, which is required by our initial condition 
that all motors started at mi = m2 = n = 0. 



A. Behavior on the filament 

For small e and small r and s, we may approximate 
(IHH) as 



1 



with J(s) 



s — ivhr + J{s) 3/(0, s) 

(83) 

which is analogous to (|42|l . This implies in real space 
that 

Pno(s) « — e-"("+-^("»/"^ (84) 

Vh 

for n > 0. Decomposing J{—s ± iO) = Ji (s) ± iJ2(s), the 
inverse Laplace transform leads to 



TTUb 



-s{t-n/vi,)-nJi{s)/vi, 



sm 



nJ2{s) 



(85) 

For large times, the variable s will be small. Since 
K{m) « (l/2)ln[16/(l - to)] for m « 1 [l|], we may 
conclude that / diverges as 

/(r,s)R,_^ln[ro(s+ir2)] with tq = -| (86) 
27r 6 16 

for small s and r. This implies in particular that 
J(s) ; 

Ji(s) ! 



31nros 

27re Inrgs 



3 [lnTos]2 + 7r2' 
and J2 s « — ^— ^- 87 

The Laplace transformed survival fraction No(s) is 
equal to Pb{0, s). It then follows from (|83|) that 



A^o(t) 



ds 



e-^V2(s) 



^ [-S + Ji(s)]2 + J|(s)' 



(88) 



For < ^ 1/e, using (|87|) and neglecting s with respect to 
Ins, we get for the number of particles on the line the 
simple result 



Noit) 



ds 3 ( 



TT 2e 



2TTet 



(89) 



which confirms the scaling No{t) ~ predicted by the 
scaling approach ^j. We can also derive the latter re- 
sult from No{s) = [3/(27re)] In(ros) using the Tauberian 
transforms summarized in appendix^ The survival frac- 
tion is shown in Fig.|51for two values of e. Again the exact 
integral (|88|) is evaluated numerically (lines) and com- 
pared to simulation data (data points). The agreement 
is good. 

Computation of the first moment A^i(s) — 
—i-§pPh{i", s)\r=o with Pb as given by leads 



to 



ds „-st [Ms)-s]J2{s) 



^ ([J,(,)_,]2 + J2(,))^ 



(90) 



For small times, this behaves again as v^t corresponding 
to the velocity Wb = l — 7 — (5— |e. For large times, we 
find 



iVi(s) 
and Ni{t) 



( 3 



V27re 
9wb 



In^(sro) 



27r2e2 



16 

3s 



(91) 



to leading order in 1/t where = 0.577215 is Euler's 
constant. This implies the average position and velocity 



nb(0 ~ — (ln-+7ij) and Vh{t) ~—^^1N^{t)vh. 
7re To 7re t 

The position of bound motors as a function of time is 
shown in Fig. EH The agreement between the analytical 
result and the simulations is again quite good. 

For the second moment, we obtain 



iV2(s) = 



2vl 



s + J(s)]3 Vh 



-Nlis) + 



[s-HJ(s)P' 
(93) 

where we have taken into account the quadratic correc- 
tion term to /(s, r = 0), and 
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ds 



e-^n3(-s + Ji(s))V2(s)- J|(s)] , 2(1-7- e) 



ds 



[(-S + Jl(s))2 + J|(s)]3 

-''*[Jl{s) - SJ,{S)J2{S) + JK-s) J2(s)] 



ds 



27re 



(94) 



The last term emerges from the singularity at s = in 
the third term of N2{s) and represents diffusion in the 
unbound state. For large times, this expression leads to 
the asymptotic relations 

3 



N2{s) « 

N2{t) « 

and 



^ \2tt( 



In'^ STo 



1 



27res' 



2vl 



2vl 



27re 

3 
27re 



i[3(lnl 

t To 

[3(ln- 



2vl 



21: e 



[(In- 

7^0 



IE) - ^ 



(95) 
27r e 



3^, (97) 



i(98) 



with To = 3/16 as in H8()|) . 
coefficient is 



The longitudinal diffusion 



H 



[In- 



To 



1e]- 



(99) 



Notice that, at typical times t ~ this is still of 

order 1/e, and, thus, much larger than the value with- 
out the unbinding mechanism. In contrast to the two- 
dimensional case discussed above, the leading term in 
three dimensions is given by the usual diffusion in the 
unbound state, but there are large logarithmic correc- 
tions of order (wb/ e )^- This can be seen in Fig.EJ where 
Arib is shown for both cases. 



1. Density profile of the bound motors 

As in the two-dimensional case, the spatial distribu- 
tion of the motors bound to the line can be derived in a 
somewhat explicit form. For n > one has 



PnMt) 



27rit)b 



with A= — [s + J(s)]-st. 



Vh 



(100) 

For large n and t, we may use the saddle point approxi- 
mation. The condition A' — yields 

s= , . (101) 

3(wbt — n) In Tos 

n this means that s is 



For very small e and v^t — n 
indeed small. 



2tt en 



3(fbi — n) In e 



(102) 



The saddle point values 



are 



A 
-A" 



271 en 
3ub Iutqs 
2t: en 



(1 
-(1 



1 



In tqS ' 
2 



21^ en 
3wbln(l/e) 
27r e n 



is^Vbhi' TqS InTos 3s^Wb lu^ £ 



(103) 
{104) 



The second derivative has a negative sign, which allows 
us to choose the contour from Sg.p. — ioo to Ss.p, + ioo, 
where Ss.p. is the saddle point value given in p01|l . The 
integration over Gaussian fluctuations yields 



Pn.oit) 



2ttA"vI 



en 



exp 



27re 



3z;b ln(l/e) 



nmb) 



As a function of n, this curve starts at 0, has a maximum 
and goes to zero at n = Vht. The apparent divergence 
near n = v\jt of the last expression is an artifact of the 
saddle point approximation. 

We can check the normalization: 



fbt/2 



dnPn,o{t). 



(106) 



It differs from the exact value by a factor of order 
■\/ln(l/ e ), which is anyhow only logarithmic. This oc- 
curs because the small n-behavior has not been treated 
properly by the saddle point approximation. 

We have compared the profiles of the bound motors 
obtained from the saddle point approximation with sim- 
ulation data, evaluating Eq. (|l()5|l with s taken from the 
numerical solution of Eq. (|1U1|I and normalizing the sad- 
dle point profile. As expected, agreement with the sim- 
ulation data is only good for large n. Therefore we have 
also taken the inverse Laplace transform of (|10()|l numer- 
ically. The result is shown in Fig. ^2 for t ~ 2000 and 
t = 10000 in comparison with simulation results. To ob- 
tain the simulated profiles, simulations were performed 
again with 5 x 10^ motors particles, most of which are, 
however, detached from the filament at the times the 
density profile was measured. While agreement is good 
for large n, there are deviations in the region around the 
maximum, which are probably due to the approximation 
used in Eq. 
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B. Behavior of the unbound motors 

Also for the motors that detached from the pinning 
Una, we can follow the steps of previous section. At large 
times, the transport properties obtained for the unbound 
motors dominate also the results that are obtained if av- 
erages are taken over all, bound and unbound, motors, 
because at large times the motors spend most of the time 
detached from the filament 



is smeared, with spread as large as the average. In d = 2 
the effect was stronger, since already the leading terms 
had different numerical prefactors, cf. Eq. Mlfl for _Db 
with Eq. ^ for D\\. 

For transversal transport the situation is much simpler 



(111) 



For large t one has small s, so s ^ J{s)- This implies 



1. Average position and dispersion of unbound motors 

From Eq. (|79|l . we now find the Fourier-Laplace trans- 
form of the distribution of unbound motors to be given 
by 



-Pub(q,r-,s) 



J{s) 



(107) 



for small r, q, and s. The term linear in r implies 



t I „ t 



Ami^it) = ^1,2 W « o - - + 1e) (112) 

and the transverse diffusion coefficient is given by 

D^-\-^,-\\^-Mt)]- (113) 

The limiting values are just those without the pinning 
line. They are reached for t ^ 1/ e • 



s s\s 



and Tiub ~ T7^(ln— +7b)- 

ZTT e To 



j{s)Y 

(108) 



The last result is just half of the value for the bound 
particles H92|) . As in the two-dimensional case, the lat- 
ter results gives also the average position of all, bound 
and unbound, motors at large times. The logarithmic 
growth of the average position and the corresponding 
time-dependent velocity confirm the predictions of scal- 
ing arguments, which give n{t) ^ \nt The time- 
dependent velocity is again given by VhNQ{t) at large 
times. 

For the second moment we get 



3s2 



N2is) 



2-7- 



2wb 

s 



(109) 

For small s or large t, this expression leads to the asymp- 
totic relations 



2. Density profile of unbound motors 

Finally, we derive an expression for the probability dis- 
tribution (or, equivalently, the density profile in the case 
of many non-interacting motors) away from the filament 
using again the saddle point approximation. From H107() 
we deduce 



dqq Jo(gm) _ J_ ^ . /^s 

(27r)2 g2 + 6s ~ 7o 27r + 6.s ~ 2?: "^'^'^ 

(114) 

for m = |m| 7^ 0, where Kq is a modified Bessel function. 
This leads to the profile 

Pn,m{t) = — I :^e^*ifo(mV6^)e-"[^+'^(^)l/"^ 

TTVh J 2,711 

(115) 

For closing the contour around the negative real axis we 
need 



ub 



s[s + J{s)Y 



1 

372 

b^[(ln_+7s)2_i], 



3 27r2e2 
t 



9?;2 
" [(In 



To 
t 

To 



IE) 



1 9^2 t 

6 + 4;rf^[l^^^+^-]' 



and 



(110) 



K,{^z) = ^H^'\-z) = -yi/f (z) = -^[Yoiz)+^Mz)]. 

(116) 



We get a result of the type 

3 f°° 
27rub Jo 
nJ2{s) 



X cos — Jo(m\/6s) - sin IllMfl Yo{mVH^n) 



where the logarithmic correction to the free diffusion co- 
efficient is half of the correction to the bound diffusion 
coefficient (1991 *). That such leading singularities have dif- 
ferent numerical prefactors in different quantities could 
have been anticipated, since the cloud of random walkers 



The saddle point approximation can be done for large 
enough m, where 



(118) 
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We then get for the distribution away from the filament 
(m^O) 



ds 



TTVy^ J 2i(6s)i/*-y/7rTO 



-A 



irv^ (6s)i/4^_2mA" 
where the saddle point values are 
1 



— Un — Vht)s + nJ(s) + mvhvQs] 

27r e 71 _ 1 
3ublnTos Inros 

2nen mvh yji/ (2s) 



(119) 



Sln^TQS Ubi — n — Vhrfiyji/ (2s) 



(120) 



and 



A" = 



2tt en 



(1 



Ss^i;^ In TqS Inros 



.,3/2 



27ren 



Ss^Ub In'^ e 



771^378 
,s3/2 



(121) 



with the saddle point value of s as given by 
2iT e nvh 



3 In tqS [wbi — n — Vhm^y3/{2s)] 



(122) 



C. Pinning line with several tracks 

As before we can assume that the line has k internal 
states, which for fc = 13 model the protofilaments on 
the microtubule. The motion on the individual tracks is 
described by a master equation analogous to Eq. J^Jl, 



(123) 



This can be analyzed as in the d = 2 case. One finds for u; ^ 



s + l-(l-7~i(5 - |e)e^'^ - 7 - ^ (^e^"^ + C(l " cost^) ' 



(124) 



The surviving fractions on the individual protofilaments 
are found by inserting r ~ 0. In the time representation 
it reads 

pj(i) = VP,f.o(<) = 7:^o(i) + ^ Ve-*J" e-i^*-f(i-— ) 

71 UJ^Q 

(125) 

The asymmetry decays as exp [— ( | e + ^ — ^ cos u;)t]. As 
in the two-dimensional situation, the asymmetry has no 
influence on the total occupation P„^o — J2j Pn,o site 
n along the microtubule, that was the subject of study 
in previous subsections. 

IV. VARIABLE STICKING PROBABILITY 

We now want to incorporate the possibility that a mo- 
tor need not bind to the line when it collides with it. We 
consider two approaches. 



A. Tubule above a two-dimensional plane 

Let us now consider a tubule located above the line 
TO = of a two-dimensional plane. We consider the fol- 
lowing jump rates: Detaching: line— ^plane i e ; attach- 
ing plane— > line ^ r/; reduced jumps in the plane away 
from position below the line: (1 — r])/A. The probability 
to jump from the tubule to the line with to = is equal 
to ^ e and the probability to jump from the line with 
TO = to the tubule is given by ^ rj. In addition, the 
motor particle jumps with probability (1 — ry)/4 from the 
line with to = to a neighboring line with to = ±1. So 
for T] = no reattachment occurs. For 77 = 1 the other 
extreme occurs: one cannot jump from the line to = 
to the rest of the plane; if initially all walkers were on 
the tubule, they will go no further than below it, but not 
wander in the plane. 

Thus we consider the master equations 
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P 

P^{t + 1) - 7P„(t) + (l-7-i5-i e)P„_i(0 + i <5P„+i(i) + i 77^^„,oW- 



(126) 
(127) 



r 



As initial condition we choose all motors on the location 
71 = of the tubule, 

P„™(0)=0 and P„(0) = ,5„,o. (128) 

The Fourier-Laplace transforms yield 

[1 + s — - COST — - cos(7]P(q,r) 

= ^ e^W-^^cosgPo(r) (129) 

[l + s- 7-(l-7-i(5-ie )e*'' - i (5e-"]P(r) 



= l + iryPo(r). 
We can integrate P(g', r) over g 

dq eP(r) — 7]CoaqPo{r) 



(130) 



^o(r) = 



27r 2 + 2s — cos r — cos q 



1 



[eP{r)~Tje-^Po{r)]. (131) 



With 



it follows that 



sinh n 
cosh /i = 2 + 2s — cos r 



Po{r) 



(132) 
(133) 



For the probability to be in the fluid no subtraction as 
in H49|) is needed. We have immediately 



P{q, r, s) = 



rje 



2e 



4s + 4 — 2 cos q — 2 cos r 



sinh fJ. ~ rj cos q 

A : — 

rje 1^ + smh fi 

4eeff y/j 

(q2 + 4s) (s + e^jj - iv^^r) 



P{r, s) 



(137) 



We can now look for the enhancement of the speed on 
the tubule; it was factor 7r/2 in Eq. H31|) . Let us assume 
that "q is small, so that there is a large time domain, 
where we may neglect it. Let us thus set 



s = r/ CT, 
Then we get 

and 



and 



1 + 2V^ 



1 + 2V^ 



cr(l + 2Vct) + i^/a 



(138) 



(139) 



.(140) 



sinh fi + rje 
Eliminating Pg now yields 

P(r, s) = (^s + (1 - 7)(1 - cosr) 

1 ri \ -1 
+ -e[cosr — —]~ivhsmr] .(134) 

2 sinh fi + rje I 

For = ^ the 77-dependent term becomes eT^ , so then 
the previous situation is recovered. For rj ~ 1 one can 
check that P{q,r,s) = Po(r, s), showing that no motors 
reach the fluid (total sticking) . 
For small parameters one has 

, , / 1 , 1 , , 1 1 - r; 

P(r, s) = ( s + - (1 - 7 - - e )r- + - e M - *«b?- 

= (s + ecff \/s - Wbr) ^ . (135) 
So the only effect is the effective detaching probability 

1 — 77 



e ^ Eeff = e ■ 



V 



(136) 



There are two domains: 

1) r < 1, > 1. Here No{t) 
Ni{t) = Vbte~^ Thus no{t) = Vht and v{t) = Wb- 
This is still a sharp profile. Although particles have de- 
tached, the remaining ones go firmly with the bare speed 

2) T ^ 1, CT ^ 1. This we discussed already in JSJl. 
The relation v{t) ~ ^ TTVhNo{t) just says that v{t) ^ Vh- 

B. Variable sticking probability 

Let us now include a variable sticking probability in 
our model: If a motor reaches the filament, it rcbinds to 
it with a probability TTad , while it is reflected from the fila- 
ment with probability 1— TTad. Such a behavior can be due 
to steric constraints, if e.g. a motor with attached bead 
diffuses close to the filament, but with the bead between 
the motor and the filament. In the long time regime the 
introduction of this additional parameter is expected to 
reduce to the probability that a motor is bound to the 
filament and thus the effective time-dependent velocity 
by a factor TTad. This has been confirmed by simulations 
for the case of random walks in open compartments j^. 
In this section we show analytically that this is indeed 
the case. 
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1. The two-dimensional case 

Let us begin with the simpler case d — 2. To include 
the sticking probability, the master equations for m — 
0, ±1 have to be modified: On the lines m = ±1 the 



rate for hopping to m = 0, i.e. to the filament is 7rad/4, 
while there is a rate (1 — TTad) /4 not to jump, equivalently 
a motor on these lines attempts to hop to the filament 
with the usual rate 1/4, but the jump is rejected with 
probability 1 — TTad . The modified master equations are 



J 



(141) 



Pn,±l{t + 1) = J Pn+l.±l + ^ Pn-l,±l + \ P«,±2 + ^Pn,0 



1 — 71'ad 



Pn,±l. 



(142) 



The equivalent of Eq. Q now contains also a term with Pi (r, s) , the Fourier-Laplace transform of the probability 
distribution Pn,i{t) along the lines with m = ±1 adjacent to the filament line: 



(1 + s - i cosq - i cos r)P(g,r,s) = 1 + ^'"^^ ^ (1 - cosg)Pi(r, s) 



,3-2e-2(5 , 1-2(5 
+ ^ 7)e' — e 



1 - e 



cos q] Ph{r,s). 



4 4 " 2 

Ph{r, s) and Pi(r, s) are related via the Fourier-Laplace transform of Eq. 1)14111 : 

Pi(r, s) = — f [1 + s - (1 - 7 - I - |)e^'- - - 7]i'b(r, s) ~ 1 

TTad V III 

Using this expression for Pi(r, s), we can proceed in the same way as above and obtain 

1 + i^^(l - e"'') 

•'I'ad ^ ' 



Pbir. s) 



[s + (1 - 7)(1 - cosr) + I cosr - ivh sinr][l + ^^;^(1 - e^^)] - fe"'" 

I 



(143) 



(144) 



(145) 



For small r and s, this leads to the asymptotic relation 

1 



Phir,s) 



(146) 



which has exactly the form of Eq. (|42|l . but with an ef- 
fective detachment rate tji^^.A- Doing the analogous cal- 
culations for the unbound motors we find 



^ub(g,r,s) 



4eV^ 



(147) 



which corresponds to Eq. H49|l again with the effective 
detachment rate e/TTad. Hence in the long time regime the 
only effect of the sticking probability TTad is a rescaling of 
the detachment rate. Thus the probability for a motor 



to be bound to the filament for large times decays as 
iVo(0 ~ 7rad/(V7i'£ t^l"^) and the average displacement 
grows as ~ (7rad/e)'\/i, i.e. both quantities are reduced 
by a factor TTad as expected from the scaling approach |^ . 

2. The three-dimensional case 



For d = 3 the calculation is completely analogous. The 
sticking probability is introduced in the same way as 
above: A motor at a neighboring site of the filament 
attempts to jump to the filament with rate 1/6 as usual, 
but the attempt is only successful with probability TTad, 
so that the motor remains at its site with probability 
(1 — 7rad)/6. Then the probability distributions P(q, r, s), 
Pi{r, s) and Ph{r, s) are related via 



^'il^, 5) = ^ ([1 + 5 - (1 - 7 - I - ^^)e'"- - V- - 7]Pb(r, s) 1 



(148) 
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and 



3 + [37 + i (5 - 67 - 4e - 3(5)e*'^ - ^ (1 - e" 



(1 - e)(cosgi +cosq2)]Phir, s) 



3 + 3s — cos r - 
(1 - 7rad)[2 - COS gi - cos gzl-Pi {r, s) 
3 + 3s — cos r — cos qi — cos q2 



cos qi — cos (72 



(149) 



From these two equations we obtain a rather complicated 
expression for Ph{r, s), which, in the hmit of small s and 
r, can be reduced to 



Ph{r, s) 



1 



s — ivhr + J(s) 

with J(s) = ' (150) 

3vrad/(r = 0, s) 

where /(r, s) is the integral (|80|l . For the unbound motors 
we find 



Pub(q, r, s) 



J(s) 



(s-zwbr + J(s))[s + f + f : 



(151) 



Both equations differ from those without the parameter 
TTad, i.e. from the case TTad = 1, only by a rescaling of 
the detachment rate, just as in the two-dimensional case 
discussed above. Therefore in this case too, the long 
time displacement and the probability to be bound to the 
filament are reduced by a factor TTad, as these quantities 
are proportional to e~^. 



were found to be in excellent agreement with results from 
Monte Carlo simulations. 

The random walks of molecular motors exhibit anoma- 
lous drift behavior. In two dimensions the average po- 
sition of both the bound and unbound motors grows as 
^ \ft at large times while in three dimensions the dis- 
placements grows only logarithmically. In addition, dif- 
fusion parallel to the filament is strongly enhanced. In 
the two-dimensional case, the diffusion coefficient has an 
anomalously high value, which is of the order (wb/e)^ 
where e denoted the small detachment probability, see 
(|41|l . In the three-dimensional case, there are large loga- 
rithmic corrections to the usual diffusion behavior, again 
of the order (wb/e)^, see 

Finally let us emphasize that similar behavior is also 
obtained for random walks in confined geometries which 
have effectively the same dimensionality These ge- 
ometries are accessible to in vitro experiments. In addi- 
tion, unbinding of motors from filaments and rebinding 
to them might also be important for the design of nan- 
otechnological devices using molecular motors as trans- 
port systems, which has been proposed by several groups 

EHHIII. 



SUMMARY AND CONCLUSIONS 



In summary, we have calculated various transport 
properties arising from the random walks of molecular 
motors. Over large length scales (^ ^^^"ca)^ molecular 
motors perform random walks which consist of alternat- 
ing sequences of directed movements along filaments and 
non-directed Brownian motion in the surrounding fluid. 
Here, we have described these walks as random walks 
on a cubic lattice and derived analytical solutions for the 
cases of a single filament in two or three dimensions using 
Fourier-Laplace transforms. We have obtained closed ex- 
pressions for the probability distributions of bound and 
unbound motors and their moments, which can be eval- 
uated numerically for all times. The asymptotic behav- 
ior at small and large times was obtained fully analyti- 
cally. In this way, we derived the fraction of bound mo- 
tors, the average position and dispersion, as well as effec- 
tive velocities and diffusion coefficients. All these results 
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APPENDIX A: TAUBERIAN THEOREMS 

The Tauberian theorems allow one to obtain the 
asymptotic behavior of a function /(f) at large times 
t from the small s behavior of its Laplace transform 
/(s) = /![/(t)], see, e.g., 0|. The following inverse 
Laplace transforms >C~^[/(s)] are used for the random 
walks of molecular motors: 
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£-M-s""lns] « -^^7 -[1 + a\nt - atlj(l + a)] 

l[l + a) 

r(i + a) 

C-^[-ln's] « l[3(lni + 7i^)'-y]. (Al) 

In these expressions, T is the Gamma function, r(2:) = t^~^e~*dt, and tp the Psi function defined by il){z) = 
d[lnr(2;)]/dz, and -je ^ 0.577215 is Euler's constant. 
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FIG. 1: Random walk of a molecular motor: The motor per- 
forms directed movement along a filament (grey rod) and un- 
binds from it after a certain walking distance. The unbound 
motor diffuses in the surrounding fluid until it rebinds to the 
filament and resumes directed motion. 
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FIG. 2: The random walks of molecular motors are modeled 
as random walks on a lattice. A line of lattice sites, which is 
indicated here as a black Une, represents a filament. Motors at 
filament sites perform a directed random walks, while motors 

at non-filament sites undergo symmetric random walks. For 
the movement in two dimensions, the jump probabilities at 
non-filament sites are 1/4 for each of the four neighbor sites, at 
filament sites, a motor steps forward with probability 1 — 7 — 
(5/2 — e/2 and backward with probability 5/2; jumps to each of 
the adjacent non-filament sites which lead to unbinding occur 
with probability e/4. The dwell probability is 7. 
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FIG. 3: Fraction of motors bound to the filament, A^o, as a 
function of time t for tlie two-dimensional case. The three 
curves correspond to e = 0.03 (circles), e = 0.05 (squares) 
and e = 0.08 (diamonds). Lines are obtained from the exact 
integral 115L the data points from Monte Carlo simulation. 
In the simulations, the other parameters have been chosen as 
7 = 0, (5 = 0.6, but the results shown here depend only on e. 
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FIG. 4: Displacement of motors as a function of time t in the 
two-dimensional case, as obtained from the exact integrals 
(lines) and Monte Carlo simulations (data points): Open cir- 
cles show the average position fib of the bound motors, dia- 
monds the average position nub of the unbound motors, and 
full circles indicate the displacement n averaged over all mo- 
tors, which interpolates between the curve for the bound mo- 
tors at small times and the curve for the unbound motors at 
large times. The parameters are e = 0.05, 7 = 0, and S — 0.6. 
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FIG. 5: Density profile P\,{n,t) — P„,o(t) of motors bound 
to tlie filament as a function of the spatial coordinate n par- 
allel to the filament for different times t. The profiles are 
normalized with respect to the probability A'o to be bound 
at the filament. The lines indicate the analytical result from 
Eq. 1451 , data points are from Monte Carlo simulations. The 
parameters are e = 0.05 and 7 = 5 = 0. 




FIG. 6; Diffusion of motors parallel to the filament in d = 2 
and d = 3: Dispersion of bound motors parallel to the fila- 
ment, An^, as a function of time t. The dotted lines indicate 
the linear behavior described by the large-time difTusion co- 
efficient Dh{t = oo) as given by Eqs. (I4H and 19911 . In two 
dimensions, this diffusion coefficient is anomalously high; in 
three dimensions, it is given by the diffusion away from the 
filament but exhibits large logarithmic corrections. The pa- 
rameters are the same as in Fig. |1] 
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FIG. 7: Diffusion of motors perpendicular to tlie filament in 
d = 2: Transverse dispersion Am^ as a function of time t. 
The parameters are the same as in Fig.^] 
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FIG. 8: Random walks probabilities for movements in three 
dimensions: At a non-filament site, a motor jumps with prob- 
ability 1/6 to each of the six neighbor sites, at a filament sites, 
a forward step has probability 1 — 7 — 5/2 — 2e/3, a backward 
step (5/2, and a jump to each of the six adjacent non-filament 
sites has probability e/6; as in two dimensions, a motor does 
not step at all with probability 7. As in Fig. |5| the line of 
black lattice sites represents the filament. The shaded ares 
with white lattice sites indicate the planes spanned by the n- 
and mi axes of the lattice perpendicular to the paper plane. 
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FIG. 9: Fraction of bound motors, A^o, as a function of time 
for the three-dimensional case. The two curves are for e = 
0.03 (circles) and e — 0.05 (diamonds). The other parameters 
are 7 = and S = 0.6. 




FIG. 10: Average position fih of the motors bound to the 
filament as a function of time t in three dimensions. The 
parameters are the same as in Fig. |5| 
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FIG. 11: Density profile Ph{n) — P„^o of motors bound to the 
filament for the three-dimensional case as a function of the 
spatial coordinate n parallel to the filament at times t = 2000 
(solid line, circles) and t = 10* (dashed line, diamonds). Lines 
are obtained by numerical evaluation of the approximate an- 
alytical expression llOOII . data points are from Monte Carlo 
simulations. The parameters are the same as in Fig. 2] 



